Deblocking of interacting particle assemblies: from pinning to jamming 
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^ A wide variety of interacting particle assemblies driven by an external force are characterized 

by a transition between a blocked and a moving phase. The origin of this deblocking transition 
can be traced back to the presence of either external quenched disorder, or of internal constraints. 
' The first case belongs to the realm of the depinning transition, which, for example, is relevant for 

, flux- lines in type II superconductors and other elastic systems moving in a random medium. The 

second case is usually included within the so-called jamming scenario observed, for instance, in many 
glassy materials as well as in plastically deforming crystals. Ifere we review some aspects of the rich 
^ , phenomenology observed in interacting particle models. In particular, we discuss front depinning, 

^ . observed when particles are injected inside a random medium from the boundary, elastic and plastic 

a' depinning in particle assemblies driven by external forces, and the rheology of systems close to the 

I ' jamming transition. We emphasize similarities and differences in these phenomena. 

> , 

C/3 , I. INTRODUCTION 

a . . . . 

a Various materials ranging from synthetic nanocrystals, magnetic colloids, charged particles in Coulomb crystals, 
I ' proteins and surfactants, or vortices in type II superconductors and in Bose-Einstein condensates, form ordered self- 
, assembled structures. This topic has attracted much interest for various fundamental and technological reasons. In 
this respect, the response of these structures to external forces of various kinds (optical, magnetic, mechanical) is of 
O particular importance Q, |^ Hf. In many cases one observes the presence of blocked phases, where the evolution of 
the system is frozen. This behavior can have different origins: when it is due to the presence of quenched disorder it 
is denoted by pinning, while when it is due to intrinsic constraints it is referred to as "jamming". In both cases, a 
sufficiently large force leads to a moving phase, through a deblocking transition. 
^ ' All these systems can be modeled by a set of interacting particles moving under the action of external forces 
sometimes in a random pinning field. For instance, superconducting vortices in thin films are pinned by vacancies 
and driven by an applied current through the Lorentz force, colloids interact by Van-der Waals or dipolar forces and 
are driven by the solvent fiow. Despite the differences in these systems, one can try to identify some common features 
• in their dynamic response. This goal has been achieved mostly through the use of numerical simulations, which have 
Cf^ ' been extensively used in the past in various contexts. Here we review the results obtained from numerical simulations 
of interacting particles, in order to provide a common framework for pinning and jamming phenomena that, despite 
their similarities, have been traditionally studied by different communities. 

The transition from a blocked to a moving phase is a central problem in the theory of non-equilibrium critical 
phenomena. Beside the large body of theoretical work devoted to the depinning transition of elastic manifolds in 
disordered media and recent theories devoted to jamming in glasses and colloids, one should also mention the theory 
of absorbing-state phase transitions 0, 0] . An absorbing state is a configuration in which the evolution of the system, 
Q ' typically a stochastic lattice model, is frozen and no longer evolves. When a suitable control parameter is changed the 
O ' system can eventually be found into an active, statistically stationary, phase. The absorbing state phase transition is 
' a second order non-equilibrium phase transition, characterized by scaling laws and critical exponents, as in ordinary 
equilibrium phase transitions. The same is true for the depinning transition, and in fact it is sometimes possible to map 
a depinning transition into and absorbing-state phase transition and vice versa. While the characterization of pinning- 
depinning as a critical phenomenon is based on a firm theoretical ground, the current theoretical understanding of 
jamming phenomena is not so advanced. Similarly, when depinning involves the generation of topological defects, one 
refers to a plastic depinning transition, but the precise meaning of the word transition is not clear at present. 

At this stage of the theoretical understanding, however, it is possible to draw an extensive common picture of these 
phenomena, in which some parts are depicted in full detail, others are less precise and some are just sketched. We 
hope that this work, by its own nature incomplete, will stimulate others to find further connections between non- 
equilibrium transitions from blocked to moving phases and possibly to formulate a complete theory encompassing all 
these phenomena. 

The paper is organized as follows: in section II we discuss the models used to describe the physics of interacting 
particle assemblies. In section III, we analyze the injection of particles in a random medium and discuss the relations 
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with front propagation and with continuum theories. Section IV is devoted to the depinning of interacting particles 
by an external force through a pinning field. We analyze the problem by increasing gradually the level of complexity, 
from the pinning of single particle to collective elastic and plastic depinning. In section V, we introduce jamming 
phenomena and discuss in detail the jamming transition observed in plastically deformed crystals, modeled by a set 
of stress driven interacting dislocations. We conclude briefly in section VI. 

II. INTERACTING PARTICLE MODELS AND THEIR PHYSICAL REALIZATION 

Several systems in nature can be modeled by a collection of interacting particles. Here we summarize the main 
features of these models and discuss some concrete examples. For simplicity, we will restrict ourselves to pairwise 
interactions between particles. In this case, identifying the particle coordinates by ri, with t = 1, ...N, we can write 
the equations of motion in general as 

i 

where m is the mass of the particles, F is a damping coefficient, and J = —W{'r) is the interparticle force derived 
from an interaction potential V . The last term represents external forces, quenched disorder, or other noise sources, 
and will be discussed in detail below. In most cases of interest, dissipation is so strong that we can safely neglect 
inertia putting m = 0. Most of the following discussion will be devoted to this overdamped limit, but occasionally we 
will discuss as well the effect of inertia. 

Depending of the particular system under study, the interparticle potential can have different forms which will 
affect the dynamics of the system. The simplest case is that of a short-range repulsive central force J{f) = fisr(|f^/^„) 
which can be characterized by its peak value K{0) and its range £,y. In a series of increasing complexity, one can 
consider non-monotonic interactions (i.e. the force can be repulsive and attractive in different ranges), long-range 
interactions (i.e. K{r) ~ r~" with < a < d -I- 2, for large r), anisotropic forces (i.e. J = fK{r)), non-central forces 
(i.e. J(r) X f 0), or different combinations of the above. 

The external force normally includes a uniform driving force F, which could be time dependent. Typical examples 
are the AC drive F{t) = Fq sin(wt) or the ramp up F{t) = ct, but the possibihties are endless. In addition, one should 
consider position dependent forces due to quenched impurities that may be present in the system. Here we will mainly 
discuss a set of Np pinning points placed randomly at position giving rise to a random force field of the type 

Fpir1^J2fMr~Rp)/U (2) 
p 

where is the range of the individual pinning forces. Normally the particular shape of the pinning potential does not 
matter as long as its range is short. One can also consider correlated disorder, such as columnar and planar defects, 
depending on the particular situation at hand. In addition, thermal effects can be included adding a a random 
uncorrelated Gaussian term /^(r, t) to the equation. 

Once the interactions of the particle systems have been specified, one should also discuss the boundary and initial 
conditions of the model. A common choice is to use periodic boundary conditions, and to place the particles in their 
zero temperature equilibrium positions (i.e. forming a crystal). Alternatively, the particles can be placed randomly 
in the system mimicking a sudden quench of the system from a disordered high temperature phase. The latter may 
give rise to an intrinsic geometrical disorder. These conditions are appropriate if one is interested in modeling the 
dynamics in the bulk of the material, without worrying about surface effects. On the other hand, in some cases 
boundary effects are at the core of the phenomena and one should then implement different initial and boundary 
conditions. This case will be discussed explicitly in the next section. Periodic boundary conditions need special care 
when interactions are long ranged, since in this case one cannot impose a cutoff to the extent of the interaction force, 
as it is often done for short range-forces. One should instead consider explicitly the interaction of the particles in a 
given finite cell with all the periodic images of the system. The infinite sum over the images can rarely be performed 
exactly and since the sum is slowly converging a simple truncation of the series gives a poor approximation and may 
induce spurious effects. To overcome this problem one can employ the Ewald summation method, originally proposed 
for Coulomb interactions, after generalizing it for the appropriate interactions involved j^. 

Finally, we would like to discuss here some physical realizations of the generic model we have discussed above. 
One of the most studied examples is the flux-line lattice in type II superconductors. In superconducting films, the 
system can be really treated as a two dimensional set of interacting particles. Conversely, for thicker superconductor 
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samples one should study a set of flexible lines. In this article, we will only consider the case of rigid vortex lines. 
The interparticle force between rigid flux-lines in the framework of the London theory is given by 

J(|r) = <f>^/(8^2^3)^^(|^-,|/^)^^ (3) 

where $o is the quantized flux carried by the vortices, Ki is a Bessel function and A is the London penetration length 
[27112^ . Notice that this is a short-range (since Ki{x) ~ exp(— x) for large x) repulsive central force, with a divergence 
of the form x^^ at short distances which is cut off by the vortex core. Instead in two dimensions, the interaction is 
long-range 

JiM) - (4) 

decaying as 1/r. In addition to interaction forces, a current j flowing in the superconductor produces a Lorentz force 
F = j X B/c acting on the vortices. 

In the case of complex fluids or soft condensed matter materials Q, which contain large polymer molecules or 
colloidal particles in a solvent whose molecules are much smaller, a generic model in which inertial terms are neglected 
often provides an effective approach towards describing such systems. The solvent is considered as a continuum 
medium, characterized by its viscosity, in which energy is dissipated as the suspended particles move through it. In 
close correspondence with their characteristic dissipative motion, the suspended particles exhibit a Brownian dynamics 
due to the random collisions with the solvent molecules. This is modeled as a random Gaussian force fj in the equations 
of motion of the form 



r,{^-v^) = j{n) + ri{n,t), (5) 

where v" is the solvent velocity that can be controlled by an externally applied flow field, and J is an elastic or 
conservative force on particle i due to deformations of long molecules or bubbles, or due to other interactions (such 
as Van der Waals, electrostatic, magnetic, and excluded-volume) among the suspended particles. The amplitude of 
the autocorrelation function (|ffip) — 2TikBT is proportional to the temperature T of the system. More sophisticated 
algorithms, in which one solves similar equations to the one represented above, have been developed to model the 
rheology of dense spherical particles accounting for hydrodynamic interactions, ellipsoidal or rod- like particles 
as well as emulsions and foams . 

Another example which is worth considering from this general point of view is a collection of dislocations in a 
thin crystalline film. Crystal dislocations are topological defects characterized by a Burgers vector b in 
the case of flux-lines, in a three dimensional crystal dislocations are deformable lines. Nevertheless, one often treats 
them in the rigid approximation, obtaining an effective two dimensional particle model, which becomes exact for 
thin samples. Dislocations produce long-range stress and strain fields in the host crystal, and experience the so- 
called Peach-Koehler force due to the overall local stress. This induces an interaction force between dislocations that 
depends on their character (edge or screw, when b is perpendicular or parallel to the corresponding dislocation axis, 
respectively but that is generally long-range, decaying as 1/r, and anisotropic. For instance, the force between 

two edge dislocations at a distance r = {x,y), and with Burgers vectors in the x direction is given by 



fib'^ x{x'^ — y^) 
27r(I-zy) {x^+y^y 



where fj, is the shear modulus, and is the Poisson ratio of the host crystal. We have only considered the x 
component since, differently from flux lines, dislocations move mainly by gliding along preferential directions, namely 
the direction of the Burgers vector. Thus, while the particle system is two dimensional, the motion is confined 
along several particular directions. This fact, together with the anisotropic character of the interaction, gives rise to 
metastable structures that act as geometric constrains for their own dynamics. In this case, the driving force is often 
an externally applied stress a which acts on the dislocations through the Peach-Koehler force F = {b ■ a) x L, where 
L is the direction of the dislocation line local tangent. 



III. GRADIENT DRIVEN DYNAMICS: FRONT INVASION 



The theoretical and experimental invest igat ion of the growth dynamics of rough interfaces has became a subject of 
great scientific interest in recent years Q, . This is clearly illustrated nowadays by the large variety of studies 
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dealing with front invasion where roughening processes take place such as flow through porous media 0, ItR ll^ 
or imbibition 1191 . fla me propagation pol l2lj. deposition processes |l4l Il5| . and flux penetration in superconducting 
materials [slIsalisLli^. 

From a macroscopic point of view, the development of modeling techniques for the description of these dynamical 
systems has been generally based on the traditional approach to transport phenomena, where the governing expressions 
are usually differential equations representing local balances of the quantity of interest (e.g., mass, momentum, flux 
of superconducting vortices, etc.) in a continuum framework. However, it is sometimes unavoidable that the process 
of front propagation takes place on a particular substrate whose structural details and/or microscopic irregularities 
cannot be properly described within a standard macroscopic formalism. On the other hand, it often happens that 
these structural features may represent key factors for the development of highly efRcient materials. This is the case, 
for example, in the field of heterogeneous catalysis, where the morphological characteristics of the catalyst pore space 
can have a dramatic influence on the accessibility of the diffusion front of reagent towards the active sites in the 
deeper parts of the porous material In the extreme situation where the porous catalyst has a microscopically 
disordered, but macroscopically non-homogeneous geometry, even a departure from the classical diffusion formalism 
may be expected "2^. For example, the so-called anomalous type of transport in self-similar (fractal) structures 
usually occurs in the form of a subdiffusive regime 25, 26] . 

The invasion of magnetic flux into a disordered type II superconductor is another problem that has recently been the 
object of intense theoretical and experimental research. As a matter of fact, the magnetization properties of type II 
superconductors have been studied for ma ny y ears, but the interest in this problem has been renewed with the discovery 
of high temperature superconductors [23, ■ The magnetization process is usually described in terms of the Bean 
model [23 and its generalizations: flux lines enter into the sample and, due to the presence of quenched disorder, 
give rise to a steady flux gradient. While the Bean model provides a consistent picture of average magnetization 
properties, such as the hysteresis loop and thermal relaxation effects (^, it does not account for local fluctuations in 
time and space. It has been recently observed that flux line dynamics is intermittent, taking place in avalanches |3lj |. 
and flux fronts are not smooth IS.'l Is^l . In particular, it has been shown that the flux front crosses over from flat 
to fractal as a function of material parameters and applied field js^ . 

A widely used modeling strategy to describe the fluctuations around the Bean state consists numerical simulations 
of interacting vortices, pinned by quenched random impurities [ssL Is^ IstI Is^. Isp . Iloj . With this approach it has 
been possible to reproduce flux profiles Is^, hysteresis |^, avalanches |36l l39l l40l| and plastic flow ^35ll40|. One of 
the aims of these studies 37j is to establish precise connections between the microscopic models and the macroscopic 
behavior, captured for instance by generalized Bean models. A different approach treats the problem at mesoscopic 
scale, describing the evolution of interacting coarse-grained units [Zll , supposed to represent the system at an 
intermediate scale. While these models give a faithful representation of several features of the problem, the connection 
with the underlying microscopic dynamics still represents a very active field of research. 



A. Invading front from an interacting particle simulation 

In this section, the main objective is to show that the gradient driven dynamics of the overdamped motion of 
interacting particles in disordered substrates can display a collective behavior that is typical of front invasion processes 
with roughened interfaces. More precisely, as we show next, the idea is to provide a basic description for these system 
in terms of a particle model and to indicate how the relevant scaling laws relating the front position and the flux 
profile to the pinning strength can be consistently extracted from numerical simulations. In an infinitely long cylinder, 
the equation of motion for an interacting particle performing an overdamped motion in a random pinning landscape 
can be written as in Eq. ^ 

Tv, = J2 M - r,) + fp[{Rp - + t), (7) 
j p 

where F is the effective viscosity, the first term on the right hand side represents the particle-particle interaction, and 
the second accounts for the interaction between particles and pinning centers. Here we consider that J is a short-range 
interaction, and /p is the force due to a pinning center, modeled as a localized trap at the position Rp, with being 
the range of the wells (typically ^p ^ A), and p = 1, Np. For example, the pinning force could be modeled in terms 
of the expression, G{x) = —fox(\x\ — 1)^, for \x\ < 1 and zero otherwise. For completeness, an uncorrelated thermal 
noise term 77 with zero mean and variance (r^) = ki,T/T is also added to the equation, but we will restrict ourselves 
to the analysis of the case T = (see Ref. [43 for the implementation of thermal noise in simulations). 

For gradient driven systems, the solution and interpretation of the MD model is essentially accomplished by the 
simultaneous numerical integration of Eq. Q for each moving particle in the system, and subsequent analysis of the flux 
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FIG. 1: In the upper panel we show a typical realization of a flux front obtained from a simulation of interacting vortices in 
a disordered landscape [i^. In the lower panel, we show the growth of the average front position as a function of time for 
different pinning force strengths. The front initially grows as t^^^ and then is pinned at a distance that scales with /o. 

front propagation for different values of the pinning strength /q. For instance, one typical simulation can involve 
up to Np = 800 000 Poisson distributed pinning centers of width = A/2 in a system of size {L^ = 800A, Ly — lOOA), 
corresponding to a density of n = 10/ A^. The number N of flux lines depends essentially on the boundary condition 
adopted in the simulation. The injection of particles into the sample is implemented through the concentration at 
t = of all particles in a small slice L' <C A, parallel to the y direction, and imposing periodic boundary conditions 
in both directions. Due to mutual repulsion in the dense zone of the slice, the particles will be pushed inside the 
material, forming a penetrating front. The front position can then be taken as the x coordinate of the most advanced 
particle in the system at different times, or one can divide the system into a grid and identify the front (see Fig.^. 



B. Non-linear diffusion 



Now we show that the front penetration due to the collective motion of interacting particles in a substrate of pinning 
centers can be described by a disordered non-linear diffusion equation jsij . The equation can be obtained performing 
a coarse-graining of Eq. lO, starting from the Fokker-Planck equation for the probability distribution of the flux line 
coordinates ....,rN,t) 



at 



(8) 



where /; is the force on the particle i given by Eq. |7J. Next, we introduce the single particle density p{f,t) = 
(^^(5^(r — fi)), where the average is done over the distribution P{ri, ....,rAr,i). The evolution of p can be directly 
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obtained from Eq. Q and is given by 

where p^'^^f^r \t) is the two-point density, whose evolution depend on the three-point density and so on. The 
simplest truncation scheme involves the approximation p*^^)(r, r ' ,t) ~ p{f,t)p{r ',t). We then coarse grain the 
equation considering length scales larger than A. This can be done expanding J in Fourier space, keeping only the 
lowest order term in and retransforming back in real space. The result reads 

J (fr'j{r~ r')p{f\t) ~ -aVp{r,t), (10) 

where a = J (Prf- J(r)/2. 

The coarse graining of the disorder term is more subtle. A straightforward elimination of short wavelength modes 
would give rise, as in the previous case, to a random force Fc{f) — ~g\In, where n is the coarse grained version of 
the microscopic density of pinning centers n(r) = X^p^^l'^" -^p) ^'^^ 9 ^ fo ■ This method can not be applied for 
short-range attractive pinning forces as the one we are investigating. In this case, short wavelength modes yield a 
macroscopic contribution to pinning that can not be neglected. Consider for instance the flow between two coarse 
grained regions: short-range microscopic pinning forces give rise to a macroscopic force that should always oppose the 
motion, while the random force derived above could in principle point in the direction of the flow. In other words, 
Fc{f) should be considered as a friction force whose direction is always opposed to the driving force Fd (in our case 
Fd — p) and whose absolute value is given by \g^n\ for \Fd\ > |gVn| and to otherwise 44j. 

Collecting all the terms, we finally obtain a disordered non-linear diffusion equation for the density of particles 

r|j = V(apVp-pi?) + fcsTVV (11) 

Any solution of Eq. (|ll|l is clearly dependent on the particular boundary condition (BC) imposed to the system. 
It is therefore important to study the effect of different BC's on the dynamics of front propagation and to show that 
the macroscopic approach based on the coarse-grained expression Eq. is compatible with the MD model. To be 
specific, we consider the situation in which the front penetration takes place in a disordered type II superconductor 
and the particles are vortices interacting according to Eq. Q ^^^^1- The following BC's correspond to different 
experimental situations 0j Ell '■ 

• (A) Constant total number of vortices. Experimentally this corresponds to an external control of the magnetic 
flux. 

• (B) Constant vortex concentration at the boundary. This case corresponds to an external control of the magnetic 
field. 

• (C) Total vortex number increasing at constant rate. This represents an external control of the flux rate. 

• (D) Boundary concentration increasing at constant rate, corresponding to a constant field rate. 

One should notice that boundary conditions can be more complicated in reality, due to complex surface barriers that 
oppose flux penetration. 

For a clean system (/o = 0) at T = 0, Eq. can be solved exactly using scaling methods [47Ll48j . In this case, 
the density profiles obey the equation 

p[x,y,t)=t-^g[x/t^) , (12) 

where x and ij} satisfy x + SV' = 1. For the BC's considered: A) x = 1/3, -0 = 1/3; B) x = 0, V' = 1/2; C) x = ^1/3, 
■0 = 2/3; and D) y = —1, = 1. These exponents are in perfect agreement with the results from MD simulations 
reported in Ref. j4a |. For instance, the data reported in Fig. ^ correspond to boundary condition (A) and correctly 
scale with '0 — 1/3. The function G{u) also depends on the boundary condition and for the case (A) is given by 
G{u) = (1 — M^)/6 for u < 1 and vanishes for u > 1. The other cases are reported in Refs. 

The presence of disorder (pinning centers) induces substantial effects on the behavior of the system that can be 
quantified in terms of the front propagation and/or the shape of the density profiles of flux lines. Depending on 
the boundary conditions, it has been observed that the front is either pinned or simply slowed down 45]. Extensive 
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numerical simulations have also been performed in Refs. to show the compatibility between the MD model 

with disorder and its coarse-grained representation, Eq. Hll|l . Moreover, by varying the parameters of this continuum 
description of the front propagation, a crossover from flat to fractal flux fronts has been detected, consistent with 
experimental observations. The value of the fractal dimension suggests that the strong disorder limit is described by 
percolation. In the weak disorder limit, we recover the analytical results derived in Refs. 



IV. EXTERNALLY DRIVEN DYNAMICS: DEPINNING AND FLOW 



In this section, we discuss the response of interacting particles to an external force in presence of quenched disorder. 
The effect of quenched disorder is first analyzed in a single particle model, which, although oversimplified, still displays 
a depinning transition. The effect of the interparticle potential on the depinning transition will be first introduced in 
the framework of the elastic theory, which breaks down for strong pinning, leading to plastic depinning. Finally, we 
briefly discuss the moving phases observed for strong driving. 



A. Single particle pinning 

A first understanding of the dynamics of driven particles in a disordered landscape can be obtained considering the 
motion of a sing le particle [HlHalES- 

We consider a collection of parabolic wells for the pinning potentials, so that 
the equation of motion for a single particle is given by 

^^ = F + Y,foiXp-x)0{\x-Xp\-^p), (13) 
p 

where Xp are the random coordinates of the pinning centers, which we assume to be non-overlapping. When the particle 
is not interacting with a pinning center, it moves with constant velocity F/T, until it enters into the attraction range 
of a pinning center. Considering this as the initial condition (i.e. x{t = 0) = Xp — ^p), we can solve the equation 

dx 

T—=F + fo{Xp-x) x<Xp + ^p, (14) 

which is given by 

x{t) =Xp + F/fp - (F/fo + ep)e-*/^. (15) 

In the limit t ^ oo the particle remains trapped as long as F < F^, = ^p/o and escapes otherwise. We can interpret 
this behavior as a depinning transition: the particle is pinned for F < Fc and moves for F > F^. 

The force- velocity diagram can be computed noting that for F > ^p/o the particle spends in each trap a time 

T^V\og{{F + F,)/{F-F,)). (16) 

The total time T to cross a system of length L with Np pinning centers is given by 

T = NpT+{L^ 2^pNp)T/F (17) 



The average velocity of the particle is thus given by 

L F 



T r(Fnp log((F + F,)/(F - F + {1 - 2^pnp)) 



(18) 



where Up — Np/L is the number of pinning centers per unit length. Close to the depinning transition Fc the velocity 
increases logarithmically v ~ l/log(Fc — F), while at higher forces the velocity is proportional to the force. This 
behavior is characteristic of the depinning transition, but in general v ^ {F — Fc)^ , where /3 is a non trivial critical 
exponent. Notice that the logarithmic behavior found above is an artifact of the discontinuous pinning force. Using 
a continuous force, one obtains instead a behavior of the type w ~ (F — Fc)^/^ |53|. We notice that the dependence 
of the exponents on the pinning potential is a peculiarity of the single particle behavior and is not observed when 
interactions come into play as we will discuss in the next section. 

Single particle models are also useful for gradient driven dynamics discussed in the previous section [4^ . In that 
case the external force F is replaced by the density gradient, which can be approximated with a simple function of 
X. For instance, in the case of boundary condition (B) (see section Fill B|l F ^ Vp ^ Pq/x, where po is the boundary 
density. This simple model yields the correct scaling behavior for the front position 461] . 
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B. Elastic depinning 



In the previous section we have only considered the interactions between moving particles and static pinning centers. 
We could expect that those results will be valid only when the interactions between particles can be neglected, as 
for instance in a very diluted system. In most situations, however, interactions between particles should be explicitly 
taken into account. In this case we will still observe a depinning transition, but its quantitative and qualitative 
features will change js^, 113, 13 Issl Is^ Is^ Issf . It is convenient to first study an interacting particle system in the 
elastic approximation, in which the pinning forces are not strong enough to break the topological properties of the 
particle system. For instance, if the particles are arranged into a crystal the external force and the disorder preserve 
the topological order and no defects are generated. 

In the pinned phase we can write the positions of the particles by their displacement vectors u{Ri) — fi — Ri, where 
fi are the coordinates in the deformed system and Ri are the equilibrium positions. The particle interaction energy 
can then be expanded in terms of the displacement field, which is assumed to be small, 

U^J2 ^(^* - '^j) - ^0 + - ^o.iRj))^^iu,3m - upiRj)), (19) 

I] a, 13 I] 

where Uq is the energy in equilibrium. One can then take the continuum limit expanding in small gradient u{R) ~ 
u{R') + {R- R') ■ Vu{r) and obtain [s^ 

V^V„.\j,'rY.E..»^J^, (20, 

where the elastic tensor Eoi^p-yS can be expressed in terms of the interparticle pair potential as 

Ea,p^s = — 7: R-yRsT: (21) 

Clearly this expansion holds as long as the sum in Eg. 1211 converges. If interactions are long range, decaying with a 
slow power law V{r) - r-(''-i+'") with cr < 2, the elastic energy can not be expanded in gradients and we have to 
live with a non-local interaction. 

Coming back to the local limit, a further simplification is obtained in Eq. H2()|l if we take into account the symmetries 
of the equilibrium system. In the case of an isotropic system, the elastic tensor has only two independent components 
and the energy reduces to 

1 [j^rh-f^Wnf^V 



V = U,^-J4-'rK^-^J +,^^J . (22, 

where K and fi are the compression and shear moduli, respectively. 

Using the elastic expression for the interparticle energy, we can rewrite the equation of motion for the particles as 

^ ^ + [K + ■ u) + F + Ur, u). (23) 

When expanding around a moving state, we obtain the same equation with an extra convective term v ■ Vw on the left 
hand side |60l l61| . These elastic equations are still impossible to solve exa ctly , but several results have been obtained 
using scaling theories and renormalization group calculations [53. l53l IsiL ISSl Isfl . ISTl ISq . In the case of long-range 
interactions the gradients are replaced by a non-local interaction kernel and the equation becomes 

^= j d''r'K^p{r-r'){up{r')-U0{r))+F + UT,u), (24) 

with Kap{r) ^ i/r'^+°' for large distances. 

A first insight on the behavior of interacting particles in a disordered media in the elastic approximation can be 
gained by collective pinning theory, originally due to Larkin |6p| . In this approach, the main energetic contributions 
entering the problem are written as a function of a lengthscale L. Minimizing the sum of elastic and disorder energies, 
one obtains the characteristic length L^. at which pinning becomes relevant. Considering a region of size L, where 
the typical displacements of the elastic medium are of the order of the range of the pinning potential, the elastic 
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energy can be estimated as E^i ^ /xL''"^^^, where for simplicity we have only considered the shear modulus /i. The 
pinning energy can be estimated in the limit of weak pinning, when the main effect comes from the fluctuations in 
the disorder: in a region of size L the typical fluctuations of the pinning energy scale as -Epm ^ —Epy'LAnp, where 
Up is the density of pins. Minimizing the total energy fiU^^^^^ — Ep^jLArhp with respect to L, we obtain the Larkin 
length 



L.-[^$-] . (25) 



^pEp 

For L < Lc the elastic medium is essentially undeformcd, while elastic deformations grow strongly for L > and 
typically the system becomes rough on large lengthscales. 

Notice that the Larkin length decreases with the disorder strength as long as d < dc = 4. The dimension dc ~ 4 
enters here for the first time, but has in fact a great importance. As can be shown by renormalization group dc is the 
upper critical dimension: for d > dc the system does not roughen and the critical behavior at the depinning transition 
is essentially mean-field like. Since the physical dimensions is usually lower than four, one is lead to think that this 
limit is without practical interest. This is not completely true, since in presence of long-range interactions dc typically 
lowers, reaching sometimes the physical dimension. In particular, for non local elastic interactions with a < 2 the 
upper critical dimension is given by dc = 2a. 

Since Lc represent the characteristic lengthscale at which pinning start to be effective, Larkin estimated the depin- 
ning threshold as the force needed to unpin a region of length Lc- This can be done comparing the energy due to the 
external force Ecxt = F^pL'^ to the pinning energy at the scale L^ obtaining 

Fc ^ Ep^pL-^l-lip = ■ (26) 

Collective pinning theory provides a qualitative description of the depinning transition and a good estimate of the 
depinning threshold, but does not allow to estimate the critical exponents and the force-velocity curve. For this one 
needs renormalization group methods, which we do not want to discuss here. We just quote some results, obtained 
for the first order in a e = 4 — c? expansion: the exponent for the velocity scaling is given by /3 = 1 — e/9, while 
the roughness exponent ruling the scaling of the typical displacements with the system size m ~ is estimated as 
C, — e/3. Other exponents follow from scaling relations. For instance, one can analyze the transient behavior of the 
average particle velocity which initially decays as a power law (u) ~ and then crosses over to a constant value, 
when F > Fc, or goes exponentially to zero otherwise (see Fig. |2Jl i51|. The exponent a can be obtained by scaling 
relation as a = P/vz^ where z is the dynamic exponent and v is the correlation length exponent |5l| |. 

Simulations have been widely used in the past to obtain numerically the value of the exponents. We should 
distinguish here the vast body of work pertaining to elastic manifolds, in which the clastic approximation is enforced 
directly in the model simulating variants of Eg. 1231 55] , from particle simulations in which the original system is studied. 
In one dimension the elastic approximation works well and the exponents measured in particle models reproduce with 
a good accuracy the results obtained for elastic manifolds, namely (3 ~ 0.25 and C, — 1.25 (see Fig. O 51, 63, 64, 6^^]. 
This is not surprising since in d = 1 topological order is trivially present. In two dimensions from the mapping with 
elastic manifold we would expect (3 ~ 0.65 and C, ~ 0.75 (5^ and at least the first exponent is reproduced by particle 
simulations 1661. 



C. Plastic depinning 

When pinning forces become stronger and/or pinning centers more diluted the topological order of the system 
typically breaks down. In this case it is not possible to describe the deformation in terms of a displacement field 
as we did in the previous section and plastic deformation should be explicitly considered. Due to these difficulties, 
a complete theoretical understanding of plastic depinning is still not available and one should rely on numerical 
simulations. A simple estimate of the conditions for the occurrence of strong pinning effects can be gained in the 
framework of collective pinning theory. When Lc ^ clq, where oq is the interparticle distance, the particles are pinned 
individually. One can elaborate these type of arguments and draw phase diagrams in terms of disorder and applied 
force [6^ . 

Here we discuss the main features of plastic depinning as observed in numerical simulations. Typically the force 
velocity curve differs drastically from the one observed when depinning is elastic, which is upward convex (see Fig.|3J). 
In plastic depinning, the curve is convex downward implying that f3 > 1. For instance Ref. 1661] finds P = 2.2 in a 
simulation of colloids. Others do not provide a value for /3 since it appears that a precise estimate of the exponent 
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FIG. 2: The decay of the velocity for A'^ = 400 interacting particles with pinning in d = 1 for different values of the force. For 
F > Fc — 0.514 the velocity crosses over to a steady value. See Ref.|5H| for details on the model. 



is not straightforward [675. The main reason is that the depinning force Fc is small and one could as well believe 
that Fc = for a large system and the pinning seen in simulations could be an artifact of finite sizes. For instance, 
the curve reported in Fig. 0| can be reasonably fit as u F^ , with f3 ~ 2.7. A small value of Fc would not change 
significantly this fit but it is difficult to discriminate between the two cases. 

The important differences between plastic and elastic depinning can be highlighted computing the velocity distri- 
bution of the moving particles. In plastic depinning the solid breaks apart: some particles are pinned by the strong 
defects while others move, typically in channels. Thus the velocity distribution develops two peaks, one around zero 
corresponding to pinned particles and one at higher velocities [g^I. Using this bimodal distribution, it is possible to 
identify the fraction of moving particles (rim), using a velocity threshold situated in the middle of the two peaks (see 
Fig. Again the data are reasonably fit by a power law F'' with 7 ~ 1.8. 

The channel structure of the dynamics has been studied thoroughly in the literature analyzing the statistical 
properties of particle trajectories |33,|43,ll3- Fig.|Slwe report an example of the trajectories in the plastic regime, 
showing the separation between pinned and moving particles. It is also possible to analyze the tearing of the lattice 
studying its topological properties (i.e. the presence of defects, such as dislocations) 67- l7lll7^ . 

To understand theoretically plastic depinning several authors have introduced simplified lattice models, which are 
possibly amenable to analytical treatment. For instance, Refs. [t^, I3| propose a model in which particles flow in a 
two dimensional rotated square lattice. When the number of particles present in a site is larger than a threshold, one 
particle is transferred to a neighboring site. The flow is directed and the outlet of each site is random but fixed. It is 
interesting to note that this model is equivalent to a sandpile model [t^ [t^ , namely a random (quenched) directed 
Manna model ItI- It^ . It should be possible to use the result obtained for the directed Manna model to solve exactly 
this model |8Ct l8lj|. The model displays a plastic depinning transition: some particles are trapped and others flow in 
channels. The number of particles belonging to the flow basin scales as n„i ~ (F — Fc)'^ , with 7 ~ 0.5. The average 
particle current (or velocity) scales instead v ^ {F — Fc)^as (3 — 1.5. It was conjectured that /3 = 1 + 7, implying 
that the velocity exponent is due to the combination of the scaling of the velocity the flowing particle Vf and the one 
of the channel size. 

A similar reasoning was proposed in Ref. (75l | for plastic depinning in a disordered XY model. In this model each 
spin depins as a single particle in a smooth pinning potential as Vs ^ {F ~ FcY^^. The number of spins depinning 
also scale with the applied force as Um ^ {F — Fc)^ , yielding (3 = 7+1/2. The exponent 7 was found to be close to d, 
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FIG. 3: The force velocity curve obtained from the steady state velocity shown in Fig. |21 with pinning in d = 1. The best fit 
yields f3 = 0.22 and = 0.514 ^ 

the dimension of the lattice, which suggests that a simple geometrical description could be possible. We notice that 
the results presented in Fig. 0] satisfy approximately the relation /3 = 1 + 7, alt houg h Fc ~ 0. 

Recently, a different approach for plastic depinning was discussed in Refs. [63, through a model of a viscoelastic 
medium. In this model depending on the parameter values one observes first order type or continuous depinning. The 
model can be approached by mean- field theory [g^ and renormalization group [63, but at present the connections 
with plastic depinning in particle systems is not clear. It is interesting to notice that a first order depinning transition, 
with hysteresis, is also expected in presence of inertia |7Clj| . Raising the force from the pinned phase towards the moving 
phase is different than decreasing it when the system is already moving. In the latter case inertia will keep the system 
in motion even beyond the "depinning threshold" , resulting in an hysteretic force velocity curve. 

To summarize, plastic depinning is characterized in general by the tearing of the elastic medium through the 
production of dislocations. Only a fraction of the particles move along channels, while the others are pinned, suggesting 
to interpret the force velocity curve by a combination of scaling of single particle velocities and channel size. However, 
the wide fluctuations of the numerical values for the exponents present in the literature and the lack of a theory does 
not allow a clear quantitative picture of the phenomenon. 

D. Moving phases, hysteresis and avalanches 

When the external driving force is large enough, particles enter into a moving phase which can be of several kinds 
depending on the type of disorder and the structure of the particle system. All these aspects are reviewed in Ref. [s^ . 
A first understanding of the dynamics of an elastic system for strong driving forces comes from a high velocity 
expansion j83j . When the system flows rapidly one can write u{r, t) ^Vt + 5u, where V is the average flow velocity. 
In a first approximation the pinning force becomes effectively a thermal like noise fp{r, u) ~ fpit, Vt), with an effective 
dynamic temperature decreasing with velocity as Td ^ ^/V. Thus one would expect that beyond a certain velocity 
the elastic system would reorder since the effect of pinning forces disappear. The discreteness and periodicity of a 
particle system modifies considerably this picture: while the displacements along the direction of motion do not feel 
the disorder due to the high velocity, this is not the case for transverse displacements. Taking into account this effect, 
Giamarchi and Le Doussal show that the system should decouple in elastic channels and predict the existence 
of threshold for transverse depinning. These features have been then observed in simulations 72]. For an extensive 
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FIG. 4: (a) The force velocity curve for — 270 interacting particles with pinning in d = 2 in the strong pinning limit in 
d = 2. The fit is a power law with /3 — 2.7. (b) The velocity distribution for the same system for different values of the 
force. Notice the bimodal structure that can be used to set a threshold and identify moving particles, (c) The average number 
of moving particles scaling with the applied force as F'^ with 7 = 1.8. 



discussion of other aspects of the dynamics we refer to Ref. |82| 

So far we have only discussed the dynamics of interacting particles in random media occurring under a constant 
applied force. When the applied force is time dependent we observe other interesting phenomena. In particular, an 
AC drive leads to hysteresis, which has been studied in detail for domain walls in ferromagnetic materials js^H^, but 
the results should carry over to a generic elastic system with disorder [s^- Hysteresis is expected in the quasistatic 
limit already at the level of a single particle model driven by an elastic spring |4J|. Close to the depinning transition, 
the hysteresis loop of the average displacement as a function of the force can be obtained analytically solving the 
equation js^ 

T{du/dt) = {Fosin{ujt) - Fcfe{F - Fc). (27) 

Finally, Ref. shows that the velocity force curve also displays hysteresis, but only in the dynamic regime (i.e. 
hysteresis is lost in the limit lu when one recovers the usual force velocity depinning curve). 

It is now well established that when the force is raised slowly towards the depinning threshold the dynamics of an 
elastic system takes place in avalanches as in self-organized criticality '7(f\. In particular, the avalanche size distribution 



13 




scales as 

P{s) = s-\f{s/s*), (28) 

where the cutoff size grows as a power law with the distance from the critical point s* ~ {Fc — F)^"^ . Here v is the 
correlation length exponent and D is the fractal dimension of the avalanche (i.e. s* ^ ) where ^ is the correlation 
length. The connection with self-organized criticality is more apparent when the elastic system is driven at "constant 
velocity" [b^. This can be achieved coupling the system elastically to a slider moving at constant velocity, or, in 
other words, replacing F by kjVt — J d'^rdtu{r^ t)). In the limit oiV ^ Q and fc ^ the system is driven precisely at 
the depinning transition [g^, |23 • For a general review on avalanches and self-organized criticality see Ref. ^3 • 



V. KINETICALLY CONSTRAINED DYNAMICS: JAMMING AND SHEAR YIELDING 



In this section, we will discuss the flow behavior of a wide class of physical systems whose dynamics is governed by 
the presence of kinematical constraints induced by both interactions and geometry. One of the main features shared 
by their dynamics is the presence of jamming, a new concept recently proposed to refer to the suppression of the 
temporal relaxation of a physical system and its corresponding ability to explore the space of configurations |89l |9C| . 
Under the action of externally applied shear stresses, these systems eventually yield and are able to flow like a viscous 
fluid. Shear yielding is thus another feature they have in common. 

As we will see in the following subsections, under stress conditions both soft materials and crystalline solids [oJl 
are susceptible to display jamming and shear yielding due to the interactions and spatial arrangement of their con- 
stituent particles in the case of soft-matter systems, or to the interactions and spatial arrangement of their topological 
defects — such as dislocations — in plastically deforming crystals. Jamming and yielding could in turn be responsible 
for the remarkably similar creep and stationary flow rheology observed experimentally in these systems, in spite of 
the big differences among the materials involved. 



A. Jamming and viscoelastic flow in soft condensed matter materials 

The phenomena of jamming and yielding control the behavior andproperties of soft materials as diverse as colloidal 
suspensions '92^ , emulsions [§3 1 foams [qJ , gels , pastes [93 , biological tissues , and other soft matter 

systems. Most of these physical systems consist of various types of soft particles closely packed into an amorphous 
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state. At such high concentrations, the individual motion of particles is drastically constrained and, as a consequence, 
soft and concentrated materials usually respond like elastic solids upon the application of low stresses. On the other 
hand, they flow like viscous fluids above the so-called yield stress value a-y, exhibiting a common rheology. The 
deformation process of amorphous polymeric materials has, for instance, received a great deal of attention !7|. The 
creep compliance curve of amorphous polymer networks above the glass transition temperature has been reported to 
closely exhibit the following behavior [ofll llfl(1 | 

Jit) = 7/f7 - JO + Ct^'^ + t/r,, (29) 

where 7 is the global strain of the material, a is the external stress, jo is the instantaneous elastic component of the 
compliance, and 77 is the viscosity. According to this behavior, the Ci^/^ term, also known as Andrade creep term, 
dominates for times much smaller than the relaxation time tc characteristic of the complex fluid, while a macroscopic 
viscous flow of the form t/r] is established after much longer times t » tc- One can deflne and measure a time 
dependent effective viscosity in the Andrade creep regime r] = a/"/ ^ t^^^, which reaches its equilibrium value at 
longer times. 

A common nonlinear rheology is also observed for higher stress values, or shear strain-rate values in the case 
of constant strain-rate experiments which are in most cases performed in this class of systems. The stress-strain 
relationships in the steady state are often described by phenomenological equations of the form a — ay + 07" 
\lQl, 102], which imply a nonlinear dependence of the stress a on the shear strain rate 7. For a Newtonian dispersion 
CTj, = and n — 1, resulting in a constant viscosity coefficient rj — CT/7. If instead ay ^ 0, the equation describes a 
Bingham fluid. When n < 1 the relation is known as the Hershel-Bulkeley law, but if Uy = 0, the equation describes 
a power-law fluid, with a shear strain rate decreasing viscosity. 

Although the understanding of yield and viscoplastic flow in these materials is difficulted by the absence of a clear 
mediating mechanism, such as the motion of dislocations in a crystal, it has been argued that the common rheology 
displayed by these general class of complex fluids might be attributed to two particular shared properties of these 
materials: structural disorder and metastability; which are characteristic features of an underlying glassy dynamics 
p,0,3] . In a few words, a glassy dynamics is associated to the slow structural relaxation which takes place when some 
parts of the system are trapped by their neighbors and hav e to surmount large energy barriers to explore further more 
favorable configurations. Molecular dynamics simulations |104L llOSl IIO6I Il07j | of glass- forming liquids and polymers 
have proved of much help in this respect. About the same time, a general "jamming" scenario was also proposed [s^ 
as a common framework to understand the mechanical behavior of a broader class of non-equilibrium physical systems 
(colloidal suspensions, supercooled liquids, foams, etc. and granular media) which, in spite of their differences, exhibit 
common properties such as slow dynamics and scaling features near the so-called jamming thresh old. 

Regarding the slow relaxation dynamics characteristic of soft glassy materials, recent experiments |l08lll09| show the 
necessity of incorporating dynamical heterogeneities for its complete description. In this respect, a new light scattering 
method, introduced in Ref. [9^, allowed to detect the intermittent dynamics of a gel formed from attractive colloids. 
The dynamics is found to be intermittent due to random rearrangements which appear to be localized in time. This 
and similar experiments strongly suggest that intermittent behavior seems to be a fundamental ingredient for the 
slow relaxation in jammed materials. 

Finally, it is interesting to point out that an empirical relation known as the "Cox-Merz rule" (llCll | quite successfully 
relates the non-linear steady rheology with the linear but frequency dependent rheological properties of polymer melts. 
In particular, the Cox-Merz rule relates the steady viscosity at a given shear rate 7, to the modulus of the dynamic 
viscosity at a frequency w = 7, i.e. 77(7) = |r7*(?x;)|. Seve ral works have a posteriori tried to theoretically justify this 
empirical relation starting from basic assumptions |lllj . Indeed, the time dependent rheology which follows from 
Eq. H29|l in the Andrade regime rj — a /-^ ^ t^/^, is dimensionally equivalent to the shear-thinning behavior 77 ~ 7~2/3 
reported in Ref. 106] for the steady nonlinear rheology of a binary Lennard- Jones mixture, and it would be consistent 
with the Cox-Merz rule for w = 1/t. The shear thinning exponent could also be related to thesubaging behavior 
observed in many soft-glassy materials 77 ~ , where t^ is the so-called waiting time and /i < 1 "112] . 

In the following section, we will show that most of the attributes discussed for the case of amorphous soft-glassy 
materials are also shared by crystalline materials like soft metals or vortex lattices deforming plastically due to the 
motion of dislocations. 



B. Dislocation jamming and viscoplastic creep deformation 

The viscoplastic deformation of crystalline solids is due to the creeping motion of dislocations driven by an externally 
applied stress [13 . Ill3l Ill4l IllSj . The study of the dynamics of these linear topological defects is a subject of 
considerable interest because of its practical importance in materials design and engineering. It is also interesting 
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FIG. 6: The strain rate relaxation for different applied stresses at T = for a system of size L = 100&. The initial density of 
edge dislocations is around 1%. The solid line is the best linear fit of the a — 0.01 curve and yields dj/dt ^ f""'®^. 



from the theoretical point of view for the many features that dislocation motion shares with other complex systems 
like, for instance, flux lines in high temperature superconductors, or some of the soft matter materials discussed above. 

At the beginning of the XX century, Andrade reported that the creep deformation of soft metals at constant 
temperature and stress grow s in time according to a power law with exponent 1/3, i.e. 7 ~ t^^^ where 7 is the 
global strain of the material |116| |. More generally, the creep deformation curve usually follows the relation 7(t) = 
7o + pt^^^ + Kt, where 70 is th e ins tantaneous plastic strain, (it^^^ is known as Andrade creep, and Kt is referred 
to as linear creep regime The same qualitative behavior has since been observed in many materials 

with rather different structures leading to the conclusion that this should be a process determined by quite general 
principles, independent of most material specific properties. Notice that the creep curve for amorphous polymer 
melts introduced in the previous subsection, fol l ows exact ly the same relation. Various arguments have been proposed 
within the dislocation literature |ll3l ill4Llll7L IllSl Ill9| | to try to explain Andrade's creep. Most of them are based 
on thermally activated processes over time (or strain) dependent barriers, however, there is still a lack of consensus 
on the basic mechanism involved in the phenomenon. 

As in the case of soft-matter systems, the plastic deformation of crystals only occurs when the externally applied 
stress overcomes a threshold value, the yield stress of the material. Above this threshold value, large-scale dislocation 
motion may take place, and a steady regime of plastic deformation is eventually established. Dislocations tend to move 
cooperatively under the action of external stress due to their mutual long-range and anisotropic elastic interactions, 
which can be attractive or repulsive. As a result of these interactions and of the resulting spatial dislocation structures 
they give rise to, self-induced constraints build up in the system and the motion of dislocations may eventually cease. 
Nevertheless, small variations of the external loading, the density, the dislocation distribution or the temperature can 
enhance dislocation motion in a discontinuous and intermittent manner ^120.] . 

In Ref . [qH , the temporal relaxation of a relatively simple dislocation dynamics model was studied through numerical 
simulation. In particular, a collection of parallel straight edge dislocations with Burgers vectors bi moving in a single 
slip system under the action of constant stress was showed to give rise to Andrade-like creep at short and intermediate 
times for a wide range of applied stresses, without invoking thermally activated processes, i.e. T = (see Fig.O. The 
strain rate, which is proportional to the density of mobile dislocations d'y/dt = biVi with Vi the velocity of each 
dislocation, decays as a power law with an exponent close to 2/3, in agreement with Andrade's observations. At larger 
times, the strain-rate was observed to cross over to a linear creep regime (i.e. to a plateau signaling a steady rate of 
deformation) whenever the applied stress is larger than a critical threshold <7ci or, otherwise, to decay exponentially 
to zero. 

These results suggested that a possible interpretation of dislocation motion and the corresponding creep laws of 
crystalline materials could also be found within the general "jamming" framework proposed to encompass the wide 
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FIG. 7: The time evolution of root-mean-square velocity of four single runs of the numerical simulations for different initial 
conditions. The applied stress value in all cases is a = 0.0075. The curves are depicted in a double logarithmic scale to 
emphasize the intermittent bursts characteristic of the creep dislocation dynamics around the yield threshold. 

variety of non-equilibrium soft and glassy materials discussed previously. When jammed, these systems are unable 
to explore phase space, but they can be unjammed by changing the stress, the density, or the temperature. The 
analogies of dislocation motion and these so-called jammed systems was further explored by considering the influences 
of dislocation multiplication, and thermal-like fluctuations on the dynamics. Dislocation multiplication favors the 
rearrangements of the system and induces a linear creep regime (floviring phase) at lower stress values, but it does not 
affect the initial power-law creep. The introduction of a finite effective temperature T had a similar effect |9l). 

The detailed analysis of the model data unveiled the dislocation microscopic dynamics in the Andrade and in the 
stationary regimes: Most dislocations are arranged into metastable structures so that the stress field they generate 
in the material is screened out on large length-scales. These structures consist of small-angle dislocation boundaries 
separating slightly misoriented crystalline blocks or far more complex dislocation arrangements. If the applied stress 
is below the yield threshold, dislocations are not able to easily explore the space of configurations to find the most 
favorable structure arrangement and they are, most of the time, trapped in metastable configurations which induce 
a jamming of the system. Around the yield threshold, a small fraction of dislocations may, however, attain a higher 
mobility and provoke several intermittent rearrangements of the whole system in the course of time. The stress field 
generated by these unsettled dislocations conserves the initial long-range character, and forces the system to continue 
evolving in time in a cooperative manner to try to reduce the internal shear stress (or minimize the elastic energy) 
by exploring further more favorable arrangements. 

In Fig. d we show the root-mean-square velocity (u^)^/^(t) = [^,- vf/N]^^^ of all the dislocations {N ^ 100 — 150) 
present in a square cell of size L = 1006 as a function of time for four single runs of the numerical simulations. Thus, 
each run represents the creep behavior of a small piece (a few nanometers big) of a macroscopic system, and starts 
from four different initial dislocation configurations obtained after letting the system relax in the absence of external 
load during a given time interval. The external shear stress applied is in all cases a = 0.0075, that is, in the vicinity 
of the critical threshold tXc. We can clearly appreciate the presence of a few intermittent burst after which (u^)^/^(t) 
slowly decreases in time. Similar burst, but either positive or negative, can also be observed in the corresponding 
strain-rate curves dj/dt (not shown). Andrade's power law creep appears as a result of the averaging process over 
many of these runs, mimicking the behavior of a much bigger system. The closer is the applied stress to the threshold 
the longer is the collective power-law motion before the system falls either in the jammed or in the moving state. 
Precisely at the critical point and for the case of an infinite system, the Andrade power-law could in principle last 
indefinitely. 
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FIG. 8: The steady strain rate value for different applied stresses in a double logarithmic scale. The solid line represents a 
cubic dependence of the form dj/dt ~ a'^ which appears to be in good correspondence with the simulation data for the higher 
stress values considered. 

C. Non-linear rheology 

Above the stress threshold, the system eventually exhibits a Unear creep regime in which the dislocations present 
in the system tend to glide in a coherent manner. The dependence of the steady strain-rate value on the external 
shear stress is shown in Fig. [S] Within the error bars, the simulation data for the higher stress values considered 
can be fit quite nicely by a cubic law dependence (see the solid line in the plot). This is an interesting result since, 
if we were to co mpa re with the nonlinear rheology characteristic of amorphous polymeric networks or other soft 
glassy materials |l06| . it would correspond to an effective shear-thinning viscosity for the dislocation ensemble which 
decreases with the strai n-rat e as ry = a{d'y /dt)^^ ~ {d-y /dt)~'^^^ . This result is in good agreement with the theoretical 
results obtained in Ref. llOfij l and compatible with the power law shear-thinning behavior 77 ~ 7^" with a — 0.5 — 1.0 
observed in many different complex fluids 0. 

To summarize, the flow of dislocation structures in crystalline solids undergoing plastic deformation shares common 
features with the time-dependent linear rheology and with the nonlinear steady rheology of soft glassy materials and, 
in particular, it seems to satisfy the empirical Cox-Merz rule. Notice, however, that the concentration of dislocations 
in the crystal needs not to be too high to warrant the presence of kinematical constraints and metastability in the 
dynamics. High concentrations could be replaced in this case by the long-range character of their mutual elastic 
interactions, that favor collective motions and rearrangements, and by their ability to form intricate extended spatial 
structures (in order to screen out the stress), that tend to glide in a coherent manner and thus can hamper their 
own relative motions driving the system to a jammed state. Further work is currently under way to try to precisely 
identify the most basic mechanism responsible for these remarkable similarities. 

VI. CONCLUSIONS 

In this paper we have discussed the collective dynamics of an assembly of interacting particles and, in particular, 
the transition from a blocked to a moving phase. Transitions of this kind are observed in different contexts and 
are due to different mechanisms. When the particles are blocked by quenched disorder, one typically refers to the 
depinning transition, which can be elastic when the medium preserves its topology through the transition, or plastic 
when topological defects, such as dislocations, are generated during the dynamics. The driving force for depinning 
can be due to an externally applied field, or could be self-generated by a density gradient, as in the case of front 
propagation. When the motion is not hindered by quenched disorder, but by intrinsic constraints one usually refers 
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to a jamming transition. 

Common features of depinning and jamming phenomena are, at the macroscopic level, the observation of a non- 
trivial steady-state force-velocity curve, scaling typically as d ~ (F — Fc)^ for F > Fc, and a transient power law 
relaxation of the velocity v ~ At the microscopic level, pinning and jamming systems are both characterized by 
a complex energy landscape, with many mctastable states. This leads to an intermittent avalanche-like response to 
external perturbations. Thus despite the different origins, pinning and jamming have several properties in common 
which could be possibly used to construct a comprehensive theory of deblocking transitions. Numerical simulations 
of interacting particles have played a major role so far to elucidate the detailed nature of some of these phenomena. 
The advancement of theoretical understanding is needed to redirect numerical simulations from a purely descriptive 
point of view to a deeper level of analysis. 
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